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This article provides an introductory review of inflation and cosmological perturbation theory. I 
begin by motivating the need for an epoch of inflation during the early stages of the radiation domi- 
nated era, and describe how inflation is typically achieved using scalar fields. Then, after an overview 
of linear cosmological perturbation theory, I derive the equations governing the perturbations, and 
outline the generation of the scalar and the tensor perturbations during inflation. I illustrate that 
slow roll inflation naturally leads to an almost scale invariant spectrum of perturbations, a predic- 
tion that seems to be in remarkable agreement with the measurements of the anisotropics of the 
cosmic microwave background. I describe the constraints from the recent observations on some of 
the more popular models of inflation. I conclude with a brief discussion on the status and certain 
prospects of the inflationary paradigm. 



I. INTRODUCTION 

A. A major drawback of the hot big bang model 

The prevailing theory about the origin and the evolu- 
tion of our universe is the so-called hot big bang model. 
The model is based on two crucial observations: the 
discovery of the expansion of the universe as character- 
ized by the Bubble's law, and the existence of an ex- 
ceedingly isotropic and a perfectly thermal Cosmic Mi- 
crowave Background (CMB) radiation. Since the energy 
density of radiation falls faster with the expansion than 
that of matter, these two observations immediately point 
to the fact that the universe has expanded from a hot and 
dense early phase, when radiation, rather than matter, 
was dominant. The transition to the more recent mat- 
ter dominated epoch occurs when the radiation density 
falls sufficiently low that the photons cease to interact 
with matter. The CMB is nothing but the relic radia- 
tion which is reaching us today from this epoch of decou- 
pling. While fairly isotropic, the CMB possesses small 
anisotropics (of about one part in 10^), and the observed 
pattern of the fluctuations in the CMB provides a direct 
snapshot of the universe at this epoch. The hot big bang 
model has been rather successful in predicting, say, the 
primordial abundances of the light elements in terms of 
a single parameter, viz. the baryon-to-photon ratio, and 
the value required to fit these observations matches the 
value that has been arrived at independently from the 
structure of the anisotropics in the CMB. Despite the 
success of the hot big bang model in explaining the re- 
sults from different observations, the model has a serious 
drawback. Under the model, the CMB photons arriving 
at us today from sufficiently widely separated directions 



*I will restrict my discussion to linear perturbation theory. Con- 
sidering higher order perturbations, as is done, say, while studying 
non-Gaussianities, though it is currently attracting a lot of atten- 
tion, is beyond the scope of this introductory review, 
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in the sky could not have interacted at the time of de- 
coupling. Nevertheless, one finds that the temperature of 
the CMB photons reaching us from any two diametrically 
opposite directions hardly differ. (For a detailed descrip- 
tion of the various successes and the few shortcomings of 
the hot big bang model, and the original references for 
the different points mentioned above, I would refer the 
reader to the following texts [1-8].) 



B. The scope and success of inflation 



Inflation — which refers to a period of accelerated ex- 
pansion during the early stages of the radiation domi- 
nated epoch — provides a satisfactory resolution to above- 
mentioned shortcoming of the hot big bang model. In 
fact, in addition to offering an elegant explanation for 
the extent of homogeneity and isotropy of the back- 
ground universe, inflation also provides an attractive 
causal mechanism to generate the inhomogeneities su- 
perimposed upon it. The inflationary epoch amplifies 
the tiny quantum fluctuations present at the beginning 
of the epoch and converts them into classical perturba- 
tions which leave their imprints as anisotropics in the 
CMB. These anisotropics in turn act as seeds for the 
formation of the large scale structures that we observe 
at the present time as galaxies and clusters of galax- 
ies. With the anisotropics in the CMB being measured 
to greater and greater precision, we have an unprece- 
dented scope to test the predictions of inflation. As I 
shall discuss, the simplest models of inflation driven by 
a single, slowly rolling scalar field, generically predict a 
nearly scale invariant spectrum of primordial perturba- 
tions, which seems to be in excellent agreement with the 
recent observations of the CMB. (For a discussion on 
these different aspects of the inflationary paradigm, in 
addition to the texts [1-8], see the following reviews [9- 



23]. 
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C. The organization of this review 

This article presents an introductory survey of infla- 
tion and linear cosmological perturbation theory, and it 
is organized as follows. In the next section, I shall de- 
scribe the main shortcoming of the hot big bang model, 
viz. the horizon problem. In Sec. Ill, I shall outline how 
inflation helps to overcome the horizon problem, and il- 
lustrate how inflation can be driven with scalar flelds. I 
shall also introduce the concept of slow roll inflation, and 
discuss the solutions to the background equations for a 
certain class of inflationary models in the slow roll ap- 
proximation. In Sec. IV, I shall present an overview of 
linear perturbation theory. I shall explain the classiflca- 
tion of the perturbations into scalars, vectors and tensors, 
and derive the equations governing these perturbations. 
I shall also discuss the behavior of the scalar perturba- 
tions in a particular limit. In Sec. V, after describing 
the generation of perturbations during inflation, I shall 
calculate the scalar and the tensor spectra that arise in 
power law as well as slow roll inflation. In Sec. VI, I shall 
discuss how such spectra compare with the recent obser- 
vations, and point out the constraints on some of the 
more popular models of inflation. Finally, in Sec. VII, I 
shall conclude with a brief discussion on the status and 
some prospects of the inflationary paradigm. 

D. Conventions and notations 

Before I get down to brass tacks, let me list the vari- 
ous conventions and notations that I shall adopt. Unless I 
mention otherwise, I shall work in (3-f l)-dimensions, and 
I shall adopt the metric signature of — , — , — ). While 
Greek indices shall denote all the spacetime coordinates, 
the Latin indices shall refer to the spatial coordinates. I 
shall set ft = c = 1, but shall display G explicitly, and de- 
flne the Planck mass to be Mp = (8 tt G) .1 shall ex- 
press the various quantities in terms of either the cosmic 
time t or the conformal time rj, as is convenient. An over- 
dot and an overprime shall denote differentiation with re- 
spect to the cosmic and the conformal time coordinates of 
the Friedmann metric that describes the expanding uni- 
verse. It is useful to note here that, for any given func- 
tion, say, /, / = {f'/a) and /= [{f"/a^) - {f'a'/a^)], 
where a is the scale factor associated with the Friedmann 
metric. Lastly, since observations indicate that the uni- 
verse has a rather small curvature, as it is usually done 
in the context of inflation, I shall work with the spatially 
flat Friedmann model. 

E. A few words on the references 

The different texts [1-8] and the variety of reviews [9- 
23] that I have already mentioned discuss the motivations 

for the inflationary scenario, the different models of in- 
flation, the formulation of cosmological perturbation the- 



ory, and also how the various models compare with the 
recent observations. In what follows, I shall typically re- 
fer to one or more of these texts or reviews. I have also 
tried to refer to the original papers. However, I should 
stress that these references are often representative and 
not necessarily exhaustive. 

II. THE HORIZON PROBLEM 

Let me begin by describing the horizon problem in the 

hot big bang model. Consider a spatially flat, smooth, 
Friedmann universe described by line element 

ds^ = dt^ - a^it) dx^ = a'^{r]) (dry^ - dx^) , (1) 

where t is the cosmic time, a{t) is the scale factor, and 
T] = J [dt/a{t)] denotes the conformal time coordinate. 
In such a background, the horizon, viz. the size of a 
causally connected region, is defined as the physical ra- 
dial distance travelled by a light ray from the big bang 
singularity at t = up to a given time t. The horizon 
can be expressed in terms of the scale factor a{t) as fol- 
lows [6, 7]: 



t 




Let me now compare the linear dimension of the forward 
and the backward light cones at the time of decoupling 
in the hot big bang model. 

If one assumes that the universe was dominated by 
non-relativistic matter from the time of decoupling t^ec 
until today to, then the physical size of the region on the 
last scattering surface from which we receive the CMB is 
given by 

to 

(^Oj^dec) = Odec J — 3 (tdec ^o) ^ , (3) 

^dec 

where Odec denotes the value of the scale factor at decou- 
pling, and I have used the fact that to 3> tdoc in arriving 
at the final expression. On the other hand, if I assume 
that the universe was radiation dominated from the big 

bang until the epoch of decoupling, then the linear di- 
mension of the horizon at decoupling turns out to be 

fdoc 

— ^ = (2 tdec) ■ (4) 

a [t) 



The ratio, say, R, of the linear dimensions of the back- 
ward and the forward light cones at decoupling then re- 
duces to [20] 
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and the last equality follows from the observational fact 
that to — 10^° years, while t^^c — 10^ years. In other 
words, the linear dimension of the backward light cone 
is about 70 times larger than the forward light cone at 
decoupling. Despite this, the CMB turns out to be ex- 
traordinarily isotropic. This is the horizon problem. 

There is another way of stating the horizon problem. 
Note that the physical wavelengths associated with the 
perturbations, say, Ap , always grow as the scale fac- 
tor, i.e. Ap (X a. In contrast, in a power law expan- 
sion of the form a{t) oc f^, the Hubble radius^, viz. 



= H ^ = {a /a) ^, goes as a^/'', so that we have 
(Ap/djj) oc a[(«~^)/«l. This implies that, when q < I— 
a condition which applies to both the radiation and 
the matter dominated epochs — the physical wavelength 
grows faster than the corresponding Hubble radius as we 
go hack in time [1]. In other words, the primordial per- 
turbations need to be correlated on scales much larger 
than the Hubble radius at sufBciently early times in or- 
der to result in the anisotropics that we observe in the 
CMB. Consequently, within the hot big bang model, any 
mechanism that is invoked to generate the primordial 
fluctuations has to be intrinsically acausal. 



III. THE INFLATIONARY SCENARIO 

Actually, apart from the horizon problem, there also 
exist a few other puzzles to which the hot big bang model 
is unable to provide a satisfactory solution. These other 
issues include, the extent to which the universe happens 
to be spatially flat today (to about one part in 10^), 
and the unacceptable density of relics that would have 
been formed when high energy symmetries were broken 
in the early universe, to name just two. Amongst all 
these issues, the horizon problem is arguably the most 
significant. Moreover, it so happens that, the inflationary 
solution to the horizon problem also aids in surmounting 
the other difficulties as well [24-29]. For these reasons, 
I shall restrict my attention to the horizon problem. In 
the following sub-sections, after illustrating how inflation 
helps in overcoming the horizon problem, I shall outline 
how inflation is typically achieved with scalar flelds. 



times. If a causal mechanism is to be responsible for the 
origin of the inhomogeneities, then, clearly, these length 
scales should be inside the Hubble scales (i.e. Ap < d^) 
in the very early stages of the universe. This will be 
possible provided we have an epoch in the early universe 
during which Ap decreases faster than the Hubble radius 
as we go hack in time, i.e. if we have [5, 23] 



At \d„ 



<0, 



which then leads to the condition that 



a > 0. 



(6) 



(7) 



In other words, the universe needs to undergo a phase 
of accelerated (inflationary) expansion during the early 
stages of the radiation dominated epoch if a physical 
mechanism is to account for the generation of the pri- 
mordial fluctuations. 

In order to illustrate these points, in Fig. 1, I have 
plotted log (length), where 'length' denotes either the 
physical wavelength of a mode or the Hubble radius, 
against log a, during inflation and the radiation domi- 
nated epochs [1]. For convenience, I have chosen to de- 
scribe inflation by the power law expansion a ocf^, with 
q > 1. In such a situation, while Ap cx a during all 
the epochs, the Hubble length behaves as a^/^ and 
during inflation and the radiation dominated epochs. Ev- 
idently, all these quantities will be described by straight 
lines in the above plot. Whereas the straight lines de- 
scribing the evolution of the physical wavelengths will 
have a unit slope, those describing the Hubble radius dur- 
ing the inflationary and the radiation dominated epochs 
will have a slope of much less than unity (say, for q ^ I) 
and two, respectively. It is then clear from the figure 
that, as we go back in time, modes that leave the Hubble 
radius during the later epochs will not be inside the Hub- 
ble radius during the early stages of the universe unless 
we have a period of inflation. 



How much inflation do we need? 



A. Bringing tlie modes inside the Hubble radius 

As I discussed above, length scales of cosmological in- 
terest today (say, 1 ^ Aq < 10^ Mpc), enter the Hubble 
radius either during the radiation or the matter domi- 
nated epochs, and are outside the Hubble radius at earlier 



^ I consider the Hubble radius rather than the horizon size since, in 
any power law expansion, the Hubble radius is equivalent to the 
horizon up to a finite multiplicative constant. Also, being a local 
quantity, the Hubble radius often proves to be more convenient 
to handle than the horizon. 



Let us now try and understand as to how much infla- 
tion we shall require in order to ensure that the forward 
lightcone from the big bang to the epoch of decoupling is 
at least as large as the backward light cone from today 
to the epoch of decoupling, thereby resolving the horizon 
problem. For simplicity, let me assume that the universe 
undergoes exponential inflation, say, from time ti to tf, 
during the early stages of the radiation dominated epoch. 
Let H be the constant Hubble scale during exponential 
inflation, and A denote the extent by which the scale fac- 
tor increases during the inflationary epoch. For A » 1, 
the dominant contribution to the size of the horizon at 
decoupling arises due to the rapid expansion during in- 
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Behaviour of Hubble 
radius in non-inflationary 
cosmology 



INFLATION 



Radiation 
domination 



log a(t) 



FIG. 1: Evolution of the physical wavelength Ap (in green) and the Hubble radius (in blue) has been plotted as a function 
of the scale factor a on a logarithmic plot during the inflationary and the radiation dominated epochs. As discussed in the text, 
in such a plot, the slope of the Hubble length is much less than unity during inflation, a feature which allows to bring around 
the modes inside the Hubble radius at a suitably early time. 



flation, and it can be evaluated to be 
M<.«.0)^a..7-|.(!ap) (lip)"^, (8) 



where I have set ti ~ H^^. In such a case, the ratio of 
the forward and the backward lightcones at the epoch of 
decouphng is given by [4] 

^ (jr) - (t^) 

and, in arriving at the final number, I have chosen H ~ 
1013 GeV. Clearly, i?, ~ 1, if A ~ 10^^. Given the scale 
factor, the amount of expansion that has occurred from 
an initial time U to a time t is usually expressed in terms 
of the number of e-folds defined as follows: 

N = JdtH = ln (10) 
ti 

Since In A = In (of/ai) ~ 60, it is often said that one re- 
quires about 60 e-folds of inflation to overcome the hori- 
zon problem^. 



C. Driving inflation with scalar fields 

If p and p denote the energy density and pressure of 
the smooth component of the matter field that is driving 
the expansion, then the Einstein's equations correspond- 
ing to the line element (1) result in the following two 
Friedmann equations for the scale factor a(t): 

= (^) f. (lla) 

where, as I had mentioned, H — {a /a) is the Hubble pa- 
rameter. It is clear from the second Friedmann equation 
that, for a to be positive, we require that (p 4- 3 p) < 0. 
Neither ordinary matter (corresponding to p = 0) nor 
radiation [which corresponds to p — (p/S)] satisfy this 
condition. In such a situation, we need to identify an- 
other form of matter to drive inflation. 

Scalar fields, which are often encountered in various 
models of high energy physics, can easily help us achieve 
the necessary condition, thereby leading to infiation. 
Consider a canonical scalar field, say, (/), described by 
the potential V{(f)). Such a scalar field is governed by the 



^ In fact, 60 e-folds is roughly an observational upper bound which 
ensures that the largest scale today is inside the Hubble radius 
during the inflationary epoch [30, 31]. Actually, the number of 
e-folds needed to resolve the horizon problem depends on the 



energy scale at which inflation takes place. For instance, if in- 
flation is assumed to occur at a rather low energy scale of, say, 
H ~ 10^ GeV, then, even 50 e-folds will suflice to surmount the 
horizon problem. 
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action 



S[cl>] = Jd' 



(12) 



with the associated stress-energy tensor being given by 



■ (13) 



The symmetries of the Friedmann background --viz. ho- 
mogeneity and isotropy — imply that the sc;alar field will 
depend only on time and, hence, the resulting stress- 
energy tensor will be diagonal. Therefore, the energy 
density p and the pressure p associated with the scalar 
field simplify to 




Y I - vi,i>) 



(14a) 
(14b) 



Moreover, from the action (12), one can arrive at the 
following equation of motion for the scalar field (p in the 
Friedmann universe: 



(15) 



where = {dV/d(f)). From the above expressions for 
p and p, one finds that the condition for infiation, viz. 
{p + 3p) < 0, reduces to 



4>^ < v{^). 



(16) 



In other words, inflation can be achieved if the potential 
energy of the scalar field dominates its kinetic energy. 

Given a V{(l)) that is motivated by a high energy model, 
the first Friedmann equation (11a) and the equation (15) 
that governs the evolution of the scalar field have to be 
consistently solved for the scale factor and the scalar field, 
with suitable initial conditions. But, using the expres- 
sions (14) for the energy density and the pressure associ- 
ated with the scalar field, the Friedmann equations (11) 
can be rewritten as 



H = - 




(17a) 
(17b) 



where, for convenience. I have set (SttG) = M~^, as I 
had defined earlier. These two equations can then be 
combined to express the scalar field and the potential 
parametrically in terms of the cosmic time t as follows [4] : 



Vit) 



V2M, Jdt )J(-h), 



(18a) 
(18b) 



If we know the scale factor a(i), these two equations allow 
us to 'reverse engineer' the potential from which such a 
scale factor can arise! Using this procedure, I shall now 
'reconstruct' the potentials of two commonly considered 
models of inflation. 

Consider the power law expansion 



a{t) = ai 



(19) 



with g > 1 corresponding to inflation, and cii being an 
arbitrary constant. On substituting this scale factor in 
equation (18a) and, upon integration, one can immedi- 
ately show that the scalar field evolves as 



(f) 



VW In 



.(39-1)9 



t 

Ml 



(20) 



where Vb is a constant of integration. The potential that 

leads to such a behavior can then be obtained using equa- 
tion (18b), and the above expressions for a{t) and </)(f). 
It is found to be [32] 



V{(j)) = Vo exp 



(21) 



In a similar fashion, it is straightforward to establish that 
the potential 



V{(1)) = (Sa^p^VMl) 



1- ^ 



(22) 

where F = ^{2 an) and k = [4 (1 - /3)//J], leads to the 
following behavior for a{t): 



a{t) — ai exp {cxt^) 



(23) 



with a > Q and < /? < 1, and ai is again some arbitrary 
constant. Since this scale factor grows faster than power 
law inflation, but slower than exponential expansion, it 
is referred to as intermediate inflation [33, 34]. 



D. Slow roll inflation 

The condition (16) that the potential energy of the 
inflaton'' dominates the kinetic energy is necessary for 
inflation to take placx;. However, inflation is guaranteed, 
if the fleld rolls slowly down the potential such that 



« T/(^). 



(24) 



Moreover, it can be ensured that the field is slowly rolling 
for a sufficiently long time (to achieve the required 60 or 
so e- folds of inflation), provided 



(25) 



It is common to refer to the scalar field that drives inflation as 
the inflaton. 
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These two conditions lead to the slow roll approxima- 
tion [35-37], which, as we shall see, allows one to con- 
struct analytical solutions, both for the background and 
the perturbations. The approximation is usually de- 
scribed in terms of what are referred to as the slow roll 
parameters. Two types of slow roll parameters — the po- 
tential slow roll parameters and the Hubble slow roll 
parameters — are often considered in the literature^. I 
shall describe these two sets of parameters below, and 
discuss the solutions for the inflaton and the scale fac- 
tor for a particular class of potentials in the slow roll 
approximation. 

1. The potential slow roll parameters 

When provided with a potential V{(j)), the potential 
slow roll approximation corresponds to rcqiiiring the fol- 
lowing two dimensionless parameters [35 37]: 

and n.-Mi{^)^ M 

where V^^ = (d^F/dt/)^), to be smaU when compared to 
unity. The two quantities and -q^ are referred to as 
the Potential Slow Roll (PSR) parameters. It is straight- 
forward to show that neglecting the kinetic energy term 
((/)^/2) in the Priedmann equation (17a) and the acceler- 
ation term (j) in the equation of motion (15) for the scalar 
field are equivalent to the smallness of these parameters. 
However, it should be emphasized that the converse is 
not true. The smallness of the PSR parameters is only a 
necessary condition, and it is not sufficient to ensure that 
these terms can indeed be ignored. The reason being that 
the PSR parameters only restrict the form of the poten- 
tial, and not the dynamics of the solutions. Even if 
and rj^ are small, there is no assurance that inflation will 
take place since the value of </> can be as large as possible. 
Therefore, in addition to the two PSR parameters being 
small, the slow roll approximation actually requires the 
additional condition that the scalar field is moving slowly 
along the attractor solution determined by the equation: 
{3H(p) = — (for a detailed discussion on this point, 
see Ref. [40]). 

In spite of this shortcoming, the PSR parameters often 
prove to be very handy. For instance, given a potential, 
they immediately allow us to determine the domains and 
the parameters of the potential that can lead to inflation. 



Though, I should add that, nowadays, it seems to be more 
common to use the following hierarchy of Hubble flow func- 
tions [38, 39]: 

with Ht being the Hubble parameter evaluated at some given 
time during inflation. 



I shall now discuss two examples to illustrate the utility 
of these parameters. Consider potentials of the form [41] 

V{cb) = {Vor), (27) 

where Vb is a constant and n > 0. Let us restrict our- 
selves to the region cp > wherein V{(j)) is positive for 
all n. It is straightforward to show that the slow roll con- 
ditions [viz. (ey,?7y) <C 1] are satisfied when (j) is much 
greater than Mp. Since inflation occurs for such large 
values of the field, these potentials are often classified 
as 'large field' models [20]. Now, consider the following 
potential which describes the pseudo-Nambu-Goldstone 
Boson 

V {<!>)= A [1 + cos Wf)], (28) 

where A and / are constants that characterize the depth 
and the width of the potential. This potential 'naturally' 
leads to inflation for values of the field that are small 
when compared to the Planck scale [42]. Hence, such 
models are usually referred to as 'small field' models. 

2. The Hubble slow roll parameters 

The Hubble Slow Roll (HSR) parameters turn out to 
be a better choice to describe the slow roll approximation 
than the PSR parameters since they do not require any 
additional conditions to be satisfied for the approxima- 
tion to be valid. The HSR parameters are called so, since 
they are defined in terms of the Hubble parameter H, 
which is treated as a function of the scalar field </> [36]. 
In such a case, we can write Eq. (17) as 

^=-{2Ml)H^, (29) 

where ff^ = (di?/d(/>). This expression can then be used 
to rewrite the first Priedmann equation (17a) as follows: 

^^-(Hf) =-(2^)' 

a relation that is referred to as the Hamilton-Jacobi for- 
mulation of inflation [36]. 

Taking H{(f)) to be the primary quantity, the dimen- 
sionless HSR parameters and ^j, are deflned as fol- 
lows [40]: 




(31) 

where H^^ = {d? H / dcf)"^) . On using Eqs. (15), (29) 
and (30), these two parameters can be written as 
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where p is the energy density associated with the scalar 
field. The following points are clear from these expres- 
sions. Firstly, ^ 1 is precisely the condition required 
for neglecting the kinetic energy term in the total energy 
of the scalar field. Secondly, the limit 5^ <C 1 corre- 
sponds to the situation wherein the acceleration term of 
the scalar field can be ignored in Eq. (15) when compared 
to the term involving the velocity. Finally, the inflation- 
ary condition a > exactly corresponds to < 1. 

It should again be emphasized that, since the smallness 
of the HSR parameters ensure that ^ is small, the HSR 
approximation implies the PSR approximation. But, the 
converse does not hold without assuming the constraint 
that the inflaton is already on the attractor [40]. 



3. Solutions in the slow roll approximation 

Note that the equation of motion of the scalar field (15) 
and the first Friedmann equation (17a) can be written in 
terms of the two HSR parameters as 



V 



3M2 



(33a) 
(33b) 



The slow roll approximation corresponds to the situa- 
tion wherein the HSR parameters and 5^ satisfy the 
following conditions: 

e„«l, (5„«1, and O [el,5l,{ej^)\ ^ e^. 

(34) 

At the leading order in the slow roll approximation, the 
equations (33) above reduce to 



V 



3M2 



and {2,H4^~ -V^. 



(35) 



Being first order differential equations, given a potential, 

these equations can be easily integrated to obtain the 
solutions to the scale factor and the scalar field in the 
slow roll limit. Let me now discuss the solutions in this 
limit to the large field models (27) that I had considered 
earlier. 

For the potential (27), when n ^ 4, in the slow roll 
limit, the solution to the scalar field is given by [8, 18] 



/,[(4-n)/2] 



it) 



A[(4-n)/2 




n (n - 4) 



Mp(t-ti), 
(36) 



whereas, when n = 4, one finds that 

4>{t) ~ 0i exp- [V(Vo/3) (4Mp) {t - ti) 



(37) 



and, in both these solutions, is a constant that denotes 
the value of the scalar field at some initial time ti. For 



all n, the scale factor can be expressed in terms of these 
solutions for the scalar field as follows: 



a{t) ~ Oi exp - 



1 



2nM2 



(38) 



with Ui being the value of the scalar factor at ti. It is 

also useful to note that, in the slow roll limit, the two 
equations (35) allow us to express the number of e-folds 
from ti to t during inflation as 



N : 



dtH 



1-) 

p / 



(39) 

where the upper limit (f> is the value of the scalar field at 
the time t. In terms of e-folds, for the large field models, 
the scalar field and the Hubble parameter are given by 



H'^{N) 



(2M» N] , 



(40a) 



^) -(2niV) 



in/2) 



(40b) 



IV. ESSENTIAL LINEAR, COSMOLOGICAL 
PERTURBATION THEORY 

Though the inflationary scenario was originally pro- 
posed to resolve the different puzzles related to the 
smooth background, it was soon realized that it also of- 
fers a simple mechanism to generate the primordial per- 
turbations. Before I turn to demonstrating how infla- 
tion produces these inhomogeneities, I shall provide an 
overview of essential cosmological perturbation theory. 
(As I had mentioned, I will restrict myself to linear per- 
turbations.) In what follows, after a discussion on the 
classification of the perturbations into scalars, vectors 
and tensors, I shall derive the equations governing these 
perturbations. Since the scalar perturbations are primar- 
ily responsible for the anisotropics in the CMB and the 
formation of structures, I shall also briefly highlight the 
evolution of these perturbations at super-Hubble scales 
during the radiation and the matter dominated epochs. 



A. Classification of the perturbations 

CMB observations indicate that the anisotropics at the 
epoch of decoupling are rather small (one part in 10^, as 
I had mentioned). If so, the amplitude of the deviations 
from homogeneity will be even smaller at earlier epochs. 
This suggests that the generation and the evolution of 
the perturbations (until structures begin to form late in 
the matter dominated epoch), can be studied using linear 
perturbation theory. 

In a Friedmann background, the metric perturbations 
can be decomposed according to their behavior under lo- 
cal rotation of the spatial coordinates on hypersurfaces of 
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constant time. This property leads to the classification of 
the perturbations as scalars, vectors and tensors [43, 44]. 
Scalar perturbations remain invariant under rotations 
(and, hence, can be said to have zero spin). As we shall 
see, these are the principal perturbations that are re- 
sponsible for the anisotropics and the inhomogeneities in 
the universe. Vector and tensor perturbations — as their 
names indicate — transform as vectors and tensors do un- 
der rotations (and, as a result, have spins of unity and 
two, respectively). The vector perturbations are gener- 
ated by rotational velocity fields and, therefore, are also 
referred to as the vorticity modes. Finally, the tensor 
perturbations describe gravitational waves, and it is im- 
portant to note that they can exist even in the absence 
of sources [45, 46]. 



1. The number of independent scalar, vector and tensor 
degrees of freedom 

Let us now carry out the exercise of counting the num- 
ber of degrees of freedom associated with these different 
types of perturbations [15]. To highlight the counting 
procedure, in this sub-section, I shall, in fact, work in ar- 
bitrary spacetime dimensions. In (Z)+ l)-dimcnsions, the 
metric tensor, being symmetric, has [{D + 1){D + 2)/2] 
degrees of freedom. However, not all of these degrees 
are independent since there exist {D + 1) degrees of free- 
dom associated with the coordinate transformations. If 
we eliminate these degrees, we are left with a total of 
[D {D + l)/2] independent degrees of freedom describing 
the metric perturbations. These degrees of freedom con- 
tain the scalar, the vector, and the tensor perturbations. 

Let Sg^j^u denote the metric perturbations in the Fried- 
mann universe. The perturbed metric can be split as 
follows: 



5gnu = {Sgoo, Sgoi, Sgij) . 



(41) 



Since it contains no running index, evidently, the pertur- 
bation Sgoo = A, say, is a scalar. As is well known, one 
can always decompose a vector into a gradient of a scalar 
and a vector that is divergence free. Hence, we can write 
Sgoi as 



Sgoi = {WiB + Si). 



(42) 



where ,B is a scalar, while Si is a vector that is divergence 
free, i.e. (ViiS*) = 0. A similar decomposition can be 
carried for the quantity Sgij by essentially repeating the 
above analysis on each of the two indices. Therefore, the 
quantity Sgij can be decomposed as 



^gij = i' 5ij + (Vj + V j Ti) 



2 



(Vi Vj + Vj Vi) 



(43) 



where ip and £ are scalar functions, J^i — like Si above — is 
a divergence free vector, and Hij is a symmetric, traceless 



and transverse tensor that satisfies the conditions HI = 
and (Vi Ti*-') = 0. Let me again count the degrees 
of freedom of the perturbed metric tensor through the 
various functions I have introduced. To describe 5gi^,y, 
we require the scalars A, B, tp and S, amounting to four 
degrees of freedom. We also need the two divergence 
free spatial vectors Si and J^i that add up to [2 {D — 1)] 
degrees. Moreover, after we impose the traceless and 
transverse conditions, the tensor Hij has 



D{D + 1) 



(D + l) 



iD + l){D-2) 



(44) 



degrees of freedom. Upon adding all these scalar, vector 
and tensor degrees, I obtain 



4 + 2(£)- 1) + 



{D + l){D-2) 



(D +1){D + 2) 



(45) 

which, as I discussed above, are the total number of de- 
grees of freedom associated with the perturbed metric in 
{D + l)-dimensions. 

In the same fashion, let me decompose the degrees of 
freedom associated with the coordinate transformations. 
In a Friedmann background, the (D+l) coordinate trans- 
formations that relate different coordinate systems that 
describe the perturbed metric can be expressed in terms 
of scalars, say, 6t and Sx, as follows: 



t^{t + 5t) and x' [x' + V' {Sx) 



(46) 



Similarly, one can construct coordinate transformations 
in terms of a divergence free vector, say, Sx"^, as 



t^t and x^ (a;* + Sx'') . 



(47) 



(It should be emphasized that these coordinate transfor- 
mations are of a particular form, and are not completely 
arbitrary. The form of these transformations is dictated 
by the fact that the difference between the coordinates 
are determined by the amplitude of the perturbations. In 
other words, it is a gauge transformation. I shall some- 
what elaborate on this point below.) There is no fur- 
ther coordinate degrees of freedom associated with the 
tensor perturbations. The two scalar quantities St and 
Sx, and the divergence free vector Sx^, constitute the 
[2 + {D — 1)] = (D -|- 1) degrees of freedom associated 
with the coordinate transformations. 

Let me now subtract the coordinate degrees of free- 
dom from the total degrees to arrive at the independent 
number of scalar, vector and tensor degrees of freedom. 
Four scalar functions, viz. {A,B,il),£), were required to 
describe the perturbed metric tensor Sgf^i,. But, there 
also exist two scalar degrees of freedom associated with 
one's choice of coordinates, i.e. St and Sx. Hence, the 
actual number of independent scalar degrees of freedom 
is (4 — 2) = 2. And, it is useful to notice that this is true 
in arbitrary spacetime dimensions. I had needed two di- 
vergence free vector functions. Si and Ti, amounting to 
a total of [2 {D — 1)] degrees of freedom, to describe the 
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perturbed metric tensor. If we subtract the (D — 1) vec- 
tor degrees of freedom corresponding to the coordinate 
transformations that can be achieved through (5x* from 
this number, one is left with [2 (£)-l)]-(D-l) = (£>-!) 
true vector degrees of freedom. As I have already pointed 
out, the tensor perturbations contain [{D + 1) {D — 2)/2] 
independent degrees of freedom. Upon adding these, I 
obtain that 



2 + (£>-!) 



{D + l){D-2) 





'D{D + iy 




2 



(48) 



which is the total number of independent degrees of free- 
dom describing the perturbed metric tensor. 

Note that, in the (3 + l)-dimensional case of our in- 
terest, there exists two each of the scalar, the vector and 
the tensor degrees of freedom. 



requirement being that the metric and the coordinates 
reduce to the standard Friedmann line element in the 
limit when the perturbations vanish. As a result, the dif- 
ference between the various coordinates have the same 
amplitude as the perturbations themselves. A particular 
choice of coordinates is called a gauge. If one changes the 
coordinate system, one would obtain another metric, i.e. 
a different gauge. The transformation from one gauge 
to another [such as Eqs. (46) and (47)] is referred to as 
a gauge transformation. There exists two approaches to 
studying the evolution of the perturbations. Either one 
constructs gauge invariant quantities, or chooses a par- 
ticular gauge and work in the specific gauge throughout. 
Being simpler and more wieldy, I shall adopt the latter 
approach. 



2. The decomposition theorem 



B. Scalar perturbations 



I have focused above on decomposing the perturbed 
metric tensor into the different types of perturbations. 
Let me now discuss the equivalent classification of the 
source of the metric perturbations, viz. the stress-energy 
tensor. Since the stress-energy tensor is a symmetric two 
tensor just as the metric tensor is, the perturbed stress- 
energy tensor can also be classified as the scalar, the vec- 
tor and the tensor components. For instance, while the 
perturbed inflaton and perfect fluids such as radiation 
and matter — are scalar sources, velocity fields with vor- 
ticity are, evidently, vector sources. Anisotropic stresses, 
when the possible scalar and vector contributions have 
been eliminated, constitute a tensor source. 

Given the perturbed metric tensor dg^i,, the corre- 
sponding Einstein tensor, say, SGf^i,, can immediately be 
calculated at the same order in the perturbations. The 
Einstein equations that relate the perturbed Einstein ten- 
sor to the perturbed stress-energy tensor, say, JT^,^, will 
then lead to the equations governing the perturbations. 
According to the decomposition theorem [7, 9], which I 
shall use without proof, at the linear order in the per- 
turbations, the scalar, the vector and the tensor pertur- 
bations decouple and, hence, they can be analyzed sepa- 
rately. In other words, each type of metric perturbation 
is aff'ected by only the same type of source. Therefore, 
the three types of perturbations can be studied indepen- 
dently of each other. 



3. Gauges 

The Friedmann line element (1) describes the expand- 
ing universe in the frame of a comoving observer. The co- 
moving observer is special since, it is with respect to such 
an observer that the universe appears homogeneous and 
isotropic. However, there is no such uniquely preferred 
frame of reference in the presence of perturbations. A 
variety of coordinate choices are possible, with the only 



I shall now turn to deriving the equations of motion 
governing the three different types of perturbations. Let 
me first consider the case of the scalar perturbations. 



1. The longitudinal gauge and the perturbed Einstein tensor 

As I mentioned above, I shall work in a particular 
gauge to describe the various perturbations. A conve- 
nient gauge to describe the scalar perturbations is the 
longitudinal gauge. If we take into account the scalar 
perturbations to the background metric (1), then, in this 
gauge, the Friedmann line element is given by 



ds^ = (1 -h 2 $) At^ - a^{t) (1 - 2 *) dx^ 



(49) 



where $ and '5 arc the two independent functions that 
describe the perturbations. (Note that this choice of 
gauge corresponds to >1 oc tA oc ^' and B = £ = Q.) At 
the linear order in the perturbations, the components of 
the perturbed Einstein tensor, viz. (5G^, corresponding 
to the line element (49) above can be evaluated to be [20] 



5Gl 



5G\ 



-QH + H^^ + (^^ V^*, (50a) 
2 Vi + (50b) 



+ (2ij- + 3fl"2) ^+ ( 4 ) V'P 



+ ( _) V'V,I?. 



(50c) 



where !? = ($-*). 
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2. Equations of motion 

The only sources of perturbations that I shall consider 

in this review will be scalar fields and perfect fluids. As 
I shall illustrate in the next section, scalar fields do not 
possess any anisotropic stress at the linear order in the 
perturbations. I shall assume that the perfect fluids that 
I consider do not contain any anisotropic stresses either. 
Under these conditions, the perturbed stress-energy ten- 
sor associated with such sources can be expressed as fol- 
lows: 

6T^ = Sp, 6T° = {ViSa) and 5T} = -5p 5], (51) 

where the quantities 5p, Sa, and Sp are the scalar quanti- 
ties that denote the perturbations in the energy density, 
the momentum flux, and the pressure, respectively. The 
first order Einstein's equations, viz. (5G(^ = (SttG) STjf, 
then lead to the required set of equations governing the 
perturbations. It is clear from the form of ^GJ above 
that, in the absence of anisotropic stresses, the corre- 
sponding Einstein equation leads to: $ = ^f. In such 
a case, the remaining three Einstein equations simplify 
to [20] 

-3 if F $) + Q = (4 TT G) 6p, (52a) 

Vi = (4 TT G) (V. da) , (52b) 

^ + 4H^+(2H + 3H'^^ ^ = {47:0) 6p. (52c) 

The first and the third of these Einstein equations can 
be combined to lead to the following differential equation 

for the Bardeen potential $ [11, 43]: 

$" + 3H (l + cl) $'-c2 

+ [2n' + {l + 3cl) n^] (47rGa2) dp'"', (53) 

where H. = (a'/a) is the conformal Hubble parameter. 
In arriving at the above equation, I have changed over to 
the conformal time coordinate, and have made use of the 
standard relation [47] 



5p= {c]^5p + 5p^'^) , 



(54) 



where = {p'/p') denotes the adiabatic speed of the per- 
turbations, and represents the non-adiabatic pres- 
sure perturbation. 



3. A conserved quantity at super-Hubble scales 

Consider the following quantity combination of the 
Bardeen potential $ and its time derivative [8, 11, 48, 49]: 



7e = $ 



2p \ + 

3n) \ p+p 



(55) 



a quantity that is referred to as the curvature perturba- 
tion^. Upon substituting this expression in Eq. (53) that 
describes the evolution of the potential and making 
use of the background equations (11), one obtains that, 
in Fourier space. 



n 



[(47rGa2) <5p^* - fc^ ^^j . (56) 



(Henceforth, the sub-scripts k shall refer to the wavenum- 

ber of the Fourier modes of the perturbations.) Note 
that, at super-Hubble scales, wherein the physical wave- 
lengths of the perturbations are much larger than the 
Hubble radius [i.e. when (k/aH) — {k/Ti) <^ 1], the 
term (c^ k^ $fe) can be neglected. If one further assumes 
that no non-adiabatic pressure perturbations are present 
(i.e. Sp'^^ = 0), say, as in the case of ideal fluids, then 
the above equation implies that TZ'/. ~ at super-Hubble 
scales. In othcir words, whc;n the perturbations arc; adia- 
batic, the curvature perturbation TZk in conserved when 
the modes are well outside the Hubble radius. 



4- Evolution of the Bardeen potential at super-Hubble scales 

Let us now make use of the conservation of curvature 
perturbation to understand how the Bardeen potential 
evolves at super-Hubble scales during the radiation and 
matter dominated eras. If I now deflne [6, 11, 18] 



where 



6 



(W2 _ n') a2 



1/2 



(57) 



(58) 



then, I find that, in Fourier space, the equation (53) that 
governs the Bardeen potential reduces to 



.2 1,2 _ I ::_ 



H 



5pl^. (59) 



In the absence of non-adiabatic pressure perturbations 
(i.e. when Jp"* = 0), the above differential equation for 
Uk simplifies to 



Uk = 0. 



(60) 



And, in the super-Hubble limit, say, as A: ^ 0, the general 
solution to this differential equation can be written as [6, 
18] 



d?7 



612 (,7) 



'- + C^{k)e{r,), (61) 



It is called so since it is proportional to the local three curvature 
on the spatial hypersurface [15, 19]. 
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where the cocfScicnts C^. and C-^ are fc-dependent con- 
stants that are determined by the initial conditions im- 
posed at early times. The corresponding Bardeen poten- 
tial $fe is then given by 



d?7 



n 



(62) 

Consider the power law expansion (19) that we had 
discussed earlier. Such an expansion can be expressed in 
terms of the conformal time coordinate as follows: 



a{r,) = {-nv) 



where 7 and H are constants given by 



7 



2q-l 



and H 



(q- l)ai 



1/9 



(63) 



(64) 



In addition to power law inflation (which corresponds 
to q > 1), the above scale factor describes the radia- 
tion and the matter dominates epochs (corresponding to 
q — (1/2) and (2/3), respectively) as well. During infla- 
tion, 7 < --2, with 7 = — 2 corresponding to exponential 
inflation (i.e. q 00). While 7 = in the case of the 
radiation dominated epoch, 7 = 1 during matter domina- 
tion. Moreover, note that, the quantity H is positive and 
—00 < 77 < during inflation, whereas, during radiation 
and matter domination, H is negative and < 77 < 00. 
It is helpful to notice that, in all these instances, — > 
corresponds to the super-Hubble limit. 

In the background (63), the quantity 6 that is defined 
in Eq. (58) is found to be 



7 + 1 



1/2 



1 



,7 + 2; \a{r]) 
so that, at super-Hubble scales, one has [cf. Eq. (62)] 



(65) 



3w + 5 



2/{3w + l) 



H[2(7+l)] 7y(2 7+3) 

where w is the following equation of state parameter: 



(66) 



«;=(p/p) = [(l-7)/3(l + 7)] 



(67) 



that is a constant in power law expansion. The first term 

in the above expression for denotes the growing mode 
(which is actually a constant), while the second term rep- 
resents the decaying mode. (Hence, the choice of sub- 
scripts G and D to the coefficient C.) Demanding finite- 
ness at very early times implies that the decaying mode 
has to be neglected, so that, at super-Hubble scales, I 
have [18] 



3(w + 1) 
3w + 5 



(68) 



This quantity vanishes when w = —I, which corresponds 
to exponential expansion driven by the cosmological con- 
stant. Therefore, it is often said that the cosmological 
constant does not induce any metric perturbations. 

Since is a constant at super-Hubble scales, in this 
limit, the curvature perturbation TZk [cf. Eq. (55)] is given 

by 



3 w 



3(m; + 1) 



'S>k^C^{k), 



(69) 



where I have made use of the expression (68) for $fe- As 
TZk is conserved and $fe is a constant at super-Hubble 
scales in power law expansion, when the modes enter 
the Hubble radius during the radiation or the matter 
dominated epochs, the Bardeen potentials at entry are 
given by [15] 



3 {w^ + 1) 



3w^ + 5 
3 {w^ + 1) 



3WM 



Ilk 



-] C,{k), (70a) 
^ ) CJfc), (70b) 



where I have made use of the fact that = (1/3) and 
u>„ = 0. These expressions also imply that, at super- 
Hubble scales, changes by a factor of (9/10) during the 
transition from the radiation to the matter dominated 
epoch [5, 15]. 

It is essentially the spcctriim of the Bardeen poten- 
tial when the modes enter the Hubble radius during the 
radiation and the matter dominated epochs that deter- 
mines the pattern of the anisotropics in the CMB and 
the structure of the universe at the largest scales that we 
observe today [5, 7, 8]. The quantity C^{k) that deter- 
mines such a spectrum has to be arrived at by solving 
the equation (59) exactly with suitable initial conditions 
imposed in the early universe. In the infiationary sce- 
nario, physically well motivated, quantum, initial condi- 
tions are imposed on the modes when they are deep inside 
the Hubble radius. As I shall illustrate in the following 
section, under these conditions, it is the background dy- 
namics during the inflationary regime that influences the 
functional form of C^{k). 



C. Vector perturbations 

As I had done in the case of the scalar perturbations, 

I shall work in a particular gauge to describe the vector 
perturbations as well. I shall choose a gauge wherein the 
Friedmann metric, when the vector perturbations have 
been included, is described by the line element [6-8] 

ds^ = di^ - a'^{t) [5ij + (Vi Fj + Vj Ft)] dx' dx^ . (71) 

(This corresponds to setting Si to zero and choosing J^i oc 

Fi.) In such a gauge, upon using the condition that Fi is 
a divergence free vector, the different components of the 
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perturbed Einstein tensor are found to be 
3 if (Vi F,- + V, .Fi) 



5G) 



(72a) 



. (72b) 



In the absence of any vector sources, according to the first 
order Einstein's equations, the non-zero components SC^ 
and 5Gj have to be equated to zero. These then imme- 
diately imply that the metric perturbation Fi vanishes 
identically. In other words, no vector perturbations are 
generated in the absence of sources with vorticity [6, 20]. 



D. Tensor perturbations 

Let us now turn to the case of the tensor perturbations. 
Upon the inclusion of these perturbations, the Friedmann 
metric can be described by the line element [5-8] 



ds^ = dt^ - a^it) {Sij + hij) dx' 



(73) 



where hij is a symmetric, transverse and traceless ten- 
sor. (Note that hij oc Tiij.) As I had discussed, the 
transverse and traceless conditions reduce the number of 
independent degrees of freedom of hij to two. These two 
degrees correspond to the two types of polarization as- 
sociated with the gravitational waves. It can be shown 
that, on imposing the transverse and the traceless condi- 
tions, the components of the perturbed Einstein tensor 
corresponding to the above line element simplify to [6] 



5Gl 



5G) 



SG'l = 0, 
2 



hij +2> H hij 



(74a) 

- 1 V^hi^ . (74b) 



In the absence of anisotropic stresses, one then arrives at 
the following differential equation describing the ampli- 
tude h of the gravitational waves [5-8] : 



h" + 2Hh' - V^ h = 0, 



(75) 



where, for later use, I have expressed the equation in 
terms of the conformal time coordinate. 



V. GENERATION OF PERTURBATIONS 
DURING INFLATION 

As I have repeatedly pointed out, the striking feature 
of inflation is the fact that it provides a natural mech- 
anism to generate the perturbations [50-53]. It is the 
quantum fluctuations associated with the inflaton that 
act as the primordial seeds for the inhomogeneities. I 
had earlier alluded to the fact that it is the spectrum of 



the Bardeen potential that determines the pattern of the 
anisotropies in the CMB and the formation of structures. 
Since the curvature perturbation is proportional to the 
Bardeen potential at super-Hubble scales [cf. Eq. (69)], 
the primary quantity of interest is the spectrum of curva- 
ture perturbations generated during inflation. The infla- 
ton being a scalar source, does not degenerate any vec- 
tor perturbations [7, 8]. However, as I have discussed, 
gravitational waves are generated even in the absence of 
sources [45, 46, 54]. The primordial gravitational waves 
are also important because they too leave their own dis- 
tinct imprints on the CMB [3, 5, 8]. In this section, I 
shall first obtain the equation of motion governing the 
curvature perturbation when the universe is dominated 
by the inflaton. I shall then quantize the curvature and 
the tensor perturbations, impose vacuum initial condi- 
tions, and evaluate the scalar and the tensor spectra in 
super-Hubble limit for the cases of power law and slow 
roll inflation. 



A. Equation of motion for the curvature 
perturbation 

As before, let cj) denote the homogeneous scalar field. 
Also, let 6(p denote the perturbation in the field. It is 
then straightforward to show that, in the metric (49), 
the components of the perturbed stress-energy tensor 
[cf. Eqs. (13) and (51)] associated with the scalar field 
can be expressed as 



5T? 



5(1) - $ + = dp (76a) 

Vi (j> Scj)^ = Vi ((5(7) , (76b) 
ST] = -(J)5^-^^^- 5^^ 5'j = -5p S'j. (76c) 



Evidently, the scalar field does not possess any 
anisotropic stress. As a result, ^ = ^ during inflation. 
On substituting the above expressions for 6p, da, and 6p 
in the first order Einstein equations (52) governing the 
scalar perturbations, one can arrive at following equation 
for the Bardeen potential: 

+ [2n' + {l + 3cl) 7f] (l-c^) V2$. (77) 

Upon comparing this equation for $ with the general 
equation (53), it is evident that the non-adiabatic pres- 
sure perturbation associated with the inflaton is given 
by 



dp"" 



1 



A-kGo? 



(78) 



In such a case, the equation (56) that describes the evo- 
lution of the curvature perturbation simplifles to 



n 



(79) 
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On differentiating this equation again with respect to 
time and, on using the background equations (17), the 
definition (55) and the Bardeen equation (77), one ob- 
tains the following equation of motion governing the 
Fourier modes of the curvature perturbation induced by 
the scalar field [8, 20]: 



21^ 

z 



0, 



(80) 



where the quantity z is given by 

z= (a^/H^ ={a(t)'/n). (81) 

It is useful to introduce the Mukhanov-Sasaki vari- 
able V that is defined as [55, 56] 



V = (Tlz). 



(82) 



The Fourier modes of this variable, say, Vk, satisfy the 
differential equation 



Vk = 0. 



(83) 



The expression on the right hand side is to be evaluated 
at super-Hubble scales [i.e. when [k/aH) = {k/Ti) <C 1] 
when the curvature perturbation approaches a constant 
value^. 

The tensor perturbations can be quantized in a simi- 
lar fashion as the curvature perturbation, and the corre- 
sponding spectrum can also be defined in the same way. 
If we write h = {u/a), then, in Fourier space, the equa- 
tion (75) describing the tensor perturbations reduces to 



Uk = 0. 



(87) 



And, it is helpful to note that this equation is essentially 
the same as the Mukhanov-Sasaki equation (83) with the 
quantity z replaced by the scale factor a. The tensor 
perturbation spectrum, say, Vrj.{k), can then be expressed 
in terms of the modes and as follows: 



V^{k) = 2 



27r2 



27r2 



\Uk\ 

a 



(88) 



B. Quantization of the perturbations and the 
definition of the power spectra 

On qTiantization, the homogeneity of the Friedmann 
background allows us to express the curvature pertur- 
bation TZ in terms of the Fourier modes TZk [satisfying 
Eq. (80)] as follows [18]: 



7t(r?,x) = 



d'^k 
(2 7r)3/2 



i k-x 



. (84) 

where the creation and the annihilation operators ak and 
al^ obey the standard commutation relations. At the lin- 
ear order in the perturbation theory that I am working 
in, the power spectrum as well as the statistical proper- 
ties of the scalar perturbations are entirely characterized 
by the two point function of the quantum field TZ. (This 
is why the perturbations generated by inflation are often 
termed as Gaussian. However, basically, this happens 
to be true since we are restricting ourselves to the lin- 
ear order in perturbation theory. I shall briefly comment 
on possible deviations from Gaussianity in the flnal sec- 
tion.) The power spectrum of the scalar perturbations, 
say, 7^s(A;), is given by the relation 

oo 

J fv,{k) ^ J ^^^^ (0|7t(^,x) 7t(^,x')|0) 



X exp-i [k- (x-x')], (85) 

where |0) is the vacuum state defined as ak|0) = V k. 
Using the decomposition (84), the scalar perturbation 
spectrum can then be obtained to be [20] 



with the expressions on the right hand sides to be evalu- 
ated in the super-Hubble limit, as in the scalar case. The 
additional factor of two in the above tensor spectrum has 
been included to take into account the two states of po- 
larization of the gravitational waves. 

The scalar and the tensor spectral indices are defined 
as [20] 



dink 



and rir. 



dlnfc 



(89) 



Notice the difference in the definition of these two quan- 
tities. Conventionally, a scale invariant scalar spectrum 
corresponds to rig = 1, while such a tensor spectrum is 
described by = 0. Finally, the tensor-to-scalar ratio r 
is defined as follows [20]: 



r{k) 



-Psik) 



(90) 



As I shall discuss below, the scalar spectral index 
and the tensor-to-scalar ratio r happen to be important 
inflationary parameters that can be constrained by the 
observations. 



® Earlier, in Subsec. IV B 3, I had illustrated that, if the non- 
adiabatic pressure perturbation 5pNA can be neglected, the cur- 
vature perturbation is conserved at super-Hubble scales. It can 
be shown that the non-adiabatic pressure perturbation associ- 
ated with the inflaton [cf. Eq. (78)] decays exponentially in the 
supcr-Hubblc limit and, hence, can be ignored [57, 58]. As a 
matter of fact, it is for this reason that the scalar perturbations 
produced by single scalar fields are usually referred to as adia- 
batic. 
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C. The Bunch-Davies initial conditions 

As I have mentioned, during inflation, the initial condi- 
tions on the perturbations arc imposed when the modes 
are well inside the Hubble scale [i.e. when {k/aH) = 
(k/n) » 1]. It is clear from Eqs. (83) and (87) that, in 
such a sub-Hubblc limit, the scalar and tensor modes Vk 
and Uk do not feel the curvature of the spacctime and, 
hence, the solutions to these modes behave in the follow- 
ing Minkowskian form: e^^*'^''-*. The assumption that 
the scalar and tensor perturbations are in the vacuum 
state them requires that Vk and are positive frequency 
modes at sub-Hubble scales, i.e. they have the asymp- 
totic form [3, 6] 



evidently, the tensor mode Uk will be described by the 
same solution as the one for v/- above. As a result, bar- 
ring an overall constant factor, the scalar and the tensor 
spectra will be of the same form. 

Upon using the series expansion of the Hankel function 
in the solution (93) and the expression (63), it can readily 
shown that the curvature perturbation TZk and the tensor 
amplitude hk approach a constant value in the super- 
Hubble limit [i.e. as {—krj) ^ 0], as expected. The scalar 
and the tensor power spectra evaluated at super-Hubble 
scales can be written as [60-63] 



2(7+2) 



(95) 



Hm {vk{v),Uk{v)) 

(fe/H)^oo 



ikri 



(91) 



It should be pointed out that the vacuum state asso- 
ciated with the modes that exhibit such a behavior is 
often referred to in the literature as the Bunch-Davies 
vacuum [59]. 



The quantities Ag and A.j, are given by 



As = 



7-h 1 



16 7r3 (7 + 2) M2 

1 \ /ir(J)| 



f \m\' 

I 2(27+1) 



7r3M2 ; V 2(27+1) 



(96a) 
(96b) 



D. The scalar and the tensor power spectra in 
power law and slow roll inflation 

In this section, my main goal will be to arrive at the 
scalar and the tensor perturbation spectra in slow roll 
inflation. However, deriving the perturbation spectra in 
the case of power law inflation proves to be highly in- 
structive in understanding the computation of the slow 
roll spectra. Therefore, before I derive the spectra in slow 
roll inflation, I shall discuss the power law case. 



1. The perturbation spectra in power law inflation 

In power law inflation, from the expressions (19) 
and (20), it is clear that 



where is the Gamma function and, as seems to be 
the convention [5, 20], I have included by hand, a factor of 
(4/M2) in A^ . Clearly, the spectral indices are constants, 
and are given by 



(n3-l)=n, = [2(7 + 2)] = - 



9- 



1 



(97) 



where I have made use of the relation (64). The resulting 
tensor-to-scalar ratio is also a constant, and is found to 
be 



16 (7 + 2) 



7 + 1 



(98) 



These results point to the fact that the scalar and tensor 
spectra turn more and more scale invariant as the scale 
factor approaches exponential inflation (i.e. as q ^ 00). 



(92) 



2. The power spectra in slow roll inflation 



Upon using the scale factor (63), the solution to the 
Mukhanov-Sasaki equation (83) that satisfies the initial 
condition (91) is found to be [60-63] 



-TT rj 



1/2 



3z[.+(i/2)](V2) Hl'^i-kr,), (93) 



where u = — [j + (1/2)], and H^^^ is the Hankel function 
of the first kind and of order u. Since, in such a power 
law case. 



Let us now evaluate the scalar and tensor spectra in 
slow roll inflation. Using equation (17b) and the expres- 
sion (32a) for the first Hubble slow roll parameter , the 
quantity z defined in equation (81) can be written as 



z = V2M^ (a^). 



(99) 



Also, note that the equations (32) defining the two Hub- 
ble slow roll parameters ej, and 6^ can be expressed in 
terms of the conformal time coordinate as follows: 



(94) 



n' 



= 1 - ( ^ ) and 5h - 



2He, 



(100) 
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Using these expressions for z and the HSR parameters, 
the term {z" / z) that appears in the Mukhanov-Sasaki 
equation (83) can be written as [13, 64, 65] 



5'. 



n 



(101) 

From the definition of above, it can also be estabUshed 
that 



(2-eJ. 



(102) 



Let us now rewrite the expression (100) above for 
as follows: 



(103) 



On integrating this expression by parts, and using the 
above definition of 5„, one obtains that 



(l-eH)W 



/ 



2 Eh (^h -'^h) 



At the loading order in the slow roll approximation [cf. 
Eq. (34)], the second term can be ignored and, at the 
same order, one can assume to be a constant. There- 
fore, we have 



(I-Eh) V 



(105) 



If we now use this expression for H in the expres- 
sions (101) and (102), then, at the leading order in the 
slow roll approximation, one gets that 



2 + 6eH - 2,5^ 

^2 



3e, 



(106a) 
(106b) 



with the slow roll parameters treated as constants. It is 
then clear from Eqs. (83) and (87) that the solutions to 
the variables and Uk will again be given in terms of 
Hankel functions as in the power law case (93), with the 
quantity v now given by 



+ 2e„ -(5, 



and i^T, 



(107) 

where the subscripts S and T refer to the scalar and 

tensor cases. 

It now remains to evaluate the two spectra in the 
super-Hubble limit. In this limit [i.e. as {—krj) 0], 

upon expanding the Hankcl function as a scries about 
the origin, the scalar and the tensor spectra can be ex- 



pressed as [20, 64] 
1 
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r(3/2). 

X 2(2-T-3) (i_e^)(2-T- 



1) 



(108b) 



where H is the Hubble parameter, and the second equal- 
ities express the asymptotic values in terms of the val- 
ues of the quantities at Hubble exit [i.e. at {—krj) = 
(1 — Cjj)^^]. (Also, I have multiplied the tensor spectrum 
by the factor of (4/M2), as I had mentioned earlier.) At 
the leading order in the slow roll approximation, the am- 
plitudes of the scalar and the tensor spectra can easily 
be read off from the above expressions. They are given 
by [20] 



IP 



Ml 
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(109a) 
(109b) 



fc=(o H) 



with the sub-scripts on the right hand side indicating 
that the quantities have to be evaluated when the modes 
cross the Hubble radius. Given a quantity, say, y, we can 
write [20] 



dlnfc 



fe=(a H) 



At 

din a 



k=(aH) 



din a 

dlnA:y^^(„jj) 
(110) 



where, in arriving at final expression, the following con- 
dition has been used: 



d In a 
dlnfc 



1, 



(111) 



k={aH) 



as H does not vary much during slow roll inflation. Using 
the expressions (109) for the power spectra, the defini- 
tions (89) of the spectral indices, and the relation (89), 

one can easily show that 

ns ~ (1 - 4 6^ + 2(5h) and n^~-{2e^). (112) 

These expressions unambiguously point to the fact the 
scalar and the tensor spectra that arise in slow roll infla- 
tion will be nearly scale invariant. The tensor-to-scalar 
ratio in the slow roll limit is found to be 



r~(16ej = -(8nj 



(113) 
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with the last equality often referred to as the consistency 
relation [13]. 



VI. COMPARISON WITH THE RECENT CMB 
OBSERVATIONS 

In this section, I shall briefly discuss as to how a power 
law primordial spectrum compares with the recent obser- 
vations of the anisotropics in the CMB. I shall also indi- 
cate how the observations constrain some of the models 
of inflation. 

The nearly scale invariant scalar power spectrum and 
the rather small tensor-to-scalar ratio that arise in slow 
roll inflation seems to be in good agreement with the 
observations of the anisotropics in the CMB. In Fig. 2, 
the angular power spectrum of the CMB temperature 
anisotropics corresponding to the concordant cosmolog- 
ical model — viz. a spatially flat, ACDM model^, and 
a nearly scale invariant primordial spectrum — has been 
plotted as a function of the multipoles. The figure also 
contains the data from the most recent observations of 
the CMB by the WMAP mission [66]. Visually, it is 

Angular Scale 
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FIG. 2: The angular power spectrum of the CMB tempera- 
ture anisotropies from the WMAP 5-year data (the black dots 
with error bars) and for the best fit concordant model, i.e. 
a spatially flat, ACDM model, with a power law primordial 
spectrum (the solid red curve). The blue band denotes the 
statistical uncertainty, known as the cosmic variance [5, 7, 8]. 
The concordant model seems to fit the data rather well. (This 
figure is from Ref. [66].) 



clear that the concordant model provides a reasonable fit 



The term ACDM model refers to the currently accepted compo- 
sition of our universe, viz. about 72% of dark energy, close to 
23% of cold (i.e. non-relativistic) dark matter, and roughly 5% of 
baryons. These numbers have been arrived at based on a variety 
of observations, including that of the CMB anisotropies [7, 8]. 



to the data. Detailed analysis of the data indicates that 
rig ~ 0.96, when the tensor contribution is completely 
ignored. Whereas, it is found that ~ 0.98 when ten- 
sors are taken into account. The data also constrains 
the tensor-to-scalar ratio to r < 0.43 at 95% Confidence 
Level (CL) [67]. 

The slow roll approximation also enables specific in- 
flationary models to be compared easily with the ob- 
servations. In order to do so, firstly, note that, in the 
slow roll limit, upon using the corresponding background 
equations (35), the scalar and tensor spectral ampli- 
tudes (109) can be expressed in terms of the potential 
V{4)) and its derivative as follows [20]: 



1 



Yl] 

^ 'P / k={aH) 



2 

3^ 
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, (114a) 



(114b) 



Secondly, during slow roll, the relations between the HSR 
and the PSR parameters can be obtained to be 



and (5„ 



('7v 



(115) 



and, as a result, for instance, ~ 1 indicates the end of 
inflation. Let us now use these expressions to arrive at 
observational constraints on the parameters of a couple 
of inflationary models. 

Observations indicate that the amplitude of the scalar 
perturbation associated with a mode that crossed the 
Hubble radius about 60 e-folds before the end of infla- 
tion is about (2 x 10~^), a constraint that is referred to 
as the COBE normalization [68]. Given a potential, the 
expressions (112)-(115) allow us to construct the scalar 
amplitude and spectral index as well as the tensor-to- 
scalar ratio in terms of the inflaton. Let me now focus 
on the large field models (27) that I had discussed ear- 
lier. In these models, inflation ends (i.e. ~ 1), when 
0c„d ^ [in/V2) Mp]. The value of the field at N e-folds 
before the end inflation can be obtained using the solu- 
tion (40a), and is given by 0„ ~ {^/[{■iN + n) n/2] M^). 
Using these expressions, it is straightforward to show 
that, for Vq = (rn^ jT) and n = 2, the COBE nor- 
malization condition leads to the constraint that m ~ 
(10-^Afp). Similarly, for = (A/4) and n = 4, the 
condition leads to A ~ 10^ 

The CMB observations also lead to useful constraints 
on the models in the n^-r plane. The various quanti- 
ties that we have obtained above allow us to express 
the scalar spectral index and the tensor-to-scalar ratio 
in terms of the number of e-folds iV counted from the end 
inflation. For the large field models, the solution (40a) 
enables us to arrive at the following expressions in the 
slow roll limit [20] 



2(n-h2) 



47V -h n 



and 
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(116) 
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I have already shown that, for the case of power law 
inflation (19) described by the potential (21), rig and r 
are constants given by Eqs. (97) and (98), respectively. 
These results for the power law case were exact. In the 
slow roll limit (i.e. when g ^ 1), they reduce to [67] 



1 



— ) and r = I — 
9/ V 9 



(117) 



These relations indicate that these models will be de- 
scribed by straight lines in the n^-r plane. In Fig. 3, 
the joint constraints on the n^-r plane from the recent 
WMAP data (and a couple of other datasets) have been 
displayed [67] . The behavior of a few inflationary models 
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FIG. 3: Joint constraints on the n^-r plane from the WMAP 
5 year data, and a couple of other data sets. The contours 
indicate the 68% and 95% CLs derived from the data. The 
top panel clearly points to the fact that the data excludes the 
n = 4 large field model at more than 95% CL for N < 60 (the 
solid line). In contrast, the n = 2 model (the dashed line) 
performs better, and falls at the boundary of the 68% region 
for N = 60. (I should again stress that the number of e-folds 
have been counted backwards from the end of inflation.) The 
bottom panel indicates the behavior of power law inflation in 
the slow roll limit (i.e. when q <^ 1). It is evident that, while 
models with q < 60 axe outside of the 95% region, models 
with q ~ 120 lie close to the boundary of the 68% region. 
(Note that p in the bottom panel is what I have called as q. 
This figure is from Ref. [67], and I would refer the reader to 
this reference for any further details.) 



have been included in the figure as well. It is evident 



from the figure that, amongst the large field models, the 
n = 2 case performs better than the n = 4 case. Also, 
in the case of the power law inflation, it is found that 
q < 60 is excluded at more than 95% CL. 



VII. STATUS AND PROSPECTS OF 
INFLATION 

I shall finally close by briefly commenting on the status 
and some prospects of the inflationary paradigm. 



A. Profusion of inflationary models 

As a broad concept, inflation can certainly be consid- 
ered as a success. However, a specific model of inflation 
that can be satisfactorily embedded in a high energy the- 
ory still eludes us. A plethora of inflationary models 
exist, quite a few of which are phenomenological in na- 
ture^. Despite the enormous amount of effort, including 
attempts at considering models that are not described 
by the canonical action, it will be fair to say that we do 
not yet have a satisfactory model. The principal hurdle 
facing the idea of inflation is to construct a model that 
is well motivated from the high energy perspective, and 
also fits the observational data well. 



B. Features in the primordial spectrum 

In this review, I had restricted myself to discussions 
on slow roll inflation, which leads to a featureless and 
nearly scale invariant primordial spectrum that seems 
to agree well with the recent observations of the CMB 
anisotropics. Though the agreement is rather good, there 
also exist a few points at the lower multipoles of the 
observed CMB angular power spectrum which lie out- 
side the cosmic variance associated with the concordant 
model. Given these observations, a handful of model in- 
dependent approaches have been constructed to recover 
the primordial power spectrum. At the smaller scales, 
all these approaches arrive at a spectrum that is nearly 
scale invariant. However, many of the approaches seem 
to unambiguously point to a sharp drop in power (along 
with a few distinct features) at the scales corresponding 
to the Hubble scale today. If future observations support 
the presence of such features in the primordial spectrum, 
then it poses an interesting challenge to build inflationary 
models that lead to the required spectrum. Conversely, 
the necessities of features will considerably restrict the 



* I would urge the reader to take a look at the figure mentioned 
in Ref. [69] , which lists the various inflationary models that have 
been considered in literature. The list is rather long, but, appar- 
ently, it is incomplete! It tells the tale. 
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class of allowed models of inflation. (For various efforts 
in this direction, see the references listed in Refs. [70, 71].) 

C. Deviations from Gaussianity 

I had earlier mentioned that the scalar perturbations 
generated during inflation are Gaussian in nature. I had 
also clarifled that this was primarily due to the fact that 
we had confined ourselves to linear perturbation theory. 
Deviations from Gaussianity can arise when takes into 
account the perturbations at the higher orders [72, 73]. 
However, the extent of the non-Gaussianity depends on 
a variety of reasons (for a recent discussion on the is- 
sue and further references, see Ref. [74]). Interestingly, 
recent re-analysis of the WMAP 5-year data seem to indi- 
cate sufficiently large non-Gaussianities (see, for instance, 
Ref. [76]). If future observations confirm such a large 
level of non-Gaussianity, then, it can result in a substan- 
tial tightening in the constraints on the inflationary mod- 
els. For example, canonical scalar field models that lead 
to a nearly scale invariant primordial spectrum contain 
only a small amount of non-Gaussianity and, hence, may 
cease to be viable [75]. However, it is known that primor- 
dial spectra with features can lead to reasonably large 
non-Gaussianities [77, 78]. Therefore, if non-Gaussianity 
indeed turns out to be large, then, either one may have 
to reconcile with the fact that the primordial spectrum 
contains features, or one possibly have to seriously con- 



sider models described by non-canonical actions, some of 
which are known to result in large Gaussianities (see, for 
instance, Refs. [79, 80]). 

Hopefully, one or both of these aspects will help us 
arrive at a satisfactory model of inflation. 
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